Molecular & Cellular Proteomics

Conventional univariate reference intervals (UniRIs) are widely used to identify abnormal biomarker values, but they evaluate each biomarker independently and do not account for coordinated deviations between biomarkers. We developed and evaluated a joint reference region (JRR) framework for plasma proteomics data using the Olink proteomics dataset generated by the UK Biobank Pharma Proteomics Project, covering approximately 3,000 plasma proteins. JRRs were estimated for selected protein pairs in a healthy reference subset, while UniRIs were estimated separately for individual proteins using the nonparametric method. Both approaches were then evaluated in ICD-defined disease subsets. Biomarker discovery revealed sparse and heterogeneous disease--protein associations, with some proteins recurring across multiple phenotypes and others showing more disease-specific patterns. The added value of JRRs varied across diseases and protein pairs. Across evaluated protein pairs, 56.5\% showed higher sensitivity under the JRR framework than the UniRI of the first protein, and 47.3\% showed higher sensitivity than the UniRI of the second protein. At the disease level, the median proportion of protein pairs with improved JRR sensitivity was 0.57. JRRs were most informative when univariate detection was limited but a subset of diseased observations was flagged only by the joint region. These findings suggest that JRRs provide a complementary approach to UniRIs by capturing abnormal joint biomarker configurations in high-dimensional proteomics data.

Parallelization of data acquisition substantially increases the throughput of mass spectrometry-based proteomics. However, parallelization also increases the density of mass spectra and consequently the overlap between ions, frustrating their analysis. To improve sequence identification and quantification from such spectra, we developed an open-source software for Joint Modeling of mass spectra (JMod). JMod models overlapping peaks as linear superpositions of their components in both MS1 and MS2 space, which permits multiplexed DIA with smaller mass offsets to increase the multiplexing capacity and thus proteomics throughput for a given plexDIA tag. This enables 9-plexDIA using 2 Da offset PSMtags, increasing throughput 9-fold while preserving quantitative accuracy and coverage depth. Furthermore, we use JMod to deconvolve simultaneous labeling by mass tags and heavy amino acids, thus increasing the throughput of metabolic pulse experiments measuring protein synthesis and degradation rates in single cells from mouse liver. By supporting enhanced decoding of highly multiplexed DIA spectra, JMod provides an open and flexible software that increases the throughput of sensitive proteomics.

Water availability is critical for plants and their microbial communities, including pathogens. The plant pathogen Phytophthora cinnamomi persists in soils with fluctuating moisture, yet cellular responses to water limitation remain poorly understood in Phytophthora and oomycetes more broadly. Although we recently characterized the proteomic response of P. cinnamomi to NaCl-induced osmotic and ionic stress, its response to PEG-mediated water limitation remains poorly understood, leaving a critical gap in our understanding of drought-relevant stress adaptation. Here, we quantified mycelial growth and profiled time-resolved proteome dynamics of P. cinnamomi during polyethylene glycol (PEG-3350)-treatment, simulating moderate water limiting conditions. Treatment with 5% PEG-3350 enhanced radial mycelial growth relative to controls, with no early growth inhibition observed. Label-free proteomics identified 1,097 protein groups, with 880 proteins shared between conditions and an asymmetric abundance profile dominated by decreasing protein abundance over time. Only a small subset of proteins increased, mainly enzymes involved in redox buffering (e.g., thioredoxin and glutaredoxin-like proteins) and mitochondrial/metabolic regulation (e.g., alternative oxidase) and mitochondrial/metabolic regulation. Hierarchical clustering revealed a potential three-phase temporal program: early translational and regulatory remodeling (1-6 HPT), sustained metabolic adjustment (6-12 HPT), and delayed engagement of redox and proteostasis functions (12-24 HPT). Network analysis demonstrated that redox-associated function was integrated throughout this adaptation, with individual clusters further specialized by cofactor preference (NADP- versus NAD-dependent enzymes) and distinct metabolic roles (malate dehydrogenase, CoA-ligase activity). This coordinated, multi-phase reorganization sustained mycelial growth despite moderate osmotic stress, indicating that P. cinnamomi employs active proteomic adaptation rather than passive stress tolerance. These findings reveal the cellular mechanisms underlying drought persistence in this invasive pathogen and suggest molecular targets for disease management under water-limited conditions.

Simultaneous quantification of monoglycosphingolipid stereoisomers is required to monitor changes in defective enzymatic pathways linked to diseases such as Gaucher Disease, Parkinson's Disease, and Krabbe Disease. Resolution of beta-glucosyl and beta-galactosyl epimers cannot be achieved by standard liquid chromatography, electrospray ionization, tandem mass spectrometry (LC-ESI-MS/MS). Separation becomes possible when field asymmetric ion mobility spectrometry (FAIMS), also known as differential mobility mass spectrometry (DMS), is added as an orthogonal separation technique to LC. FAIMS/DMS separates epimeric ion clusters in a high versus low electric field (separation voltage, SV) then redirects the target epimeric ions to the mass spectrometer through the application of a direct current (compensation voltage, CoV). Resolving SVs and CoVs must be manually determined for each lipid. Manual derivation is a labour-intensive process that requires pure synthetic standards, limiting the number of stereoisomers a user can include in an assay. To address this problem, we introduce here intelligent DMS (iDMS). iDMS is an in silico supervised neural network model that learns the ion mobility relationships between SV and CoV and the monoglycosphingolipid structural features of sugar headgroup, N-acyl chain length, and N-acyl degree of unsaturation. iDMS predicts the SV and CoV combinations capable of resolving any stereoisomer pair from a training dataset of composed of measured signal intensities across a range of SVs and CoVs of 12 lipids. This machine learning alternative to manual DMS optimization promises to accelerate the deployment of multiple-reaction-monitoring mode (MRM) RPLC-ESI-DMS-MS/MS assays for the routine and rapid quantification of biologically relevant monoglycosphingolipid stereoisomers.

Spatial omics technologies resolve molecular expression and spatial architecture at single-cell resolution, but profiling whole slides remains costly. In practice, only a few regions of interest (ROIs) are profiled, leaving the rest of the tissue unmeasured. S2-omics was the first framework to unify ROI selection with out-of-ROI prediction, but it operates on superpixels rather than individual cells and predicts discrete cell types rather than continuous molecular profiles. Superpixel-based representations do not explicitly preserve cell boundaries, while categorical cell-type labels cannot quantify molecular expression within cells. Here we present RECON, a two-stage framework that performs ROI inference and whole-slide molecular reconstruction at single-cell resolution, predicting both continuous molecular profiles and discrete cell-type labels. In the first stage, RECON extracts morphological and microenvironmental features from individual cells to identify a representative ROI for spatially resolved single-cell molecular profiling. In the second stage, RECON trains deep learning models on molecular measurements acquired within the selected ROI and reconstructs transcriptomic or proteomic profiles for all remaining cells on the slide. Benchmarked against pathologist annotations, RECONs ROI selection outperforms the superpixel-based S2-omics approaches (IoU: 0.75 versus 0.64). For transcriptomics, refining the modeling unit from superpixels to single cells improves per-gene Pearson correlation by 22%. For proteomics, RECON surpasses the current state-of-the-art method, ROSIE, across all 16 markers, with a median per-cell Pearson correlation of 0.91 versus 0.84. Moreover, RECON delineates tumour boundaries and regions with distinct immune-cell densities, and highlights candidate tertiary lymphoid structures. Together, these results demonstrate that RECON enables informative ROI selection and whole-slide molecular reconstruction at single-cell resolution for both spatial transcriptomics and spatial proteomics.

Proteins that distinguish quiescent cells from other non-proliferative states and actively regulate their return to proliferation remain poorly understood. Here, we combined quantitative proteomics with functional image-based screening to identify regulators of the quiescence-to-proliferation transition. Amongst the functional quiescence signature proteins we identified, we focussed on integrin 11 (ITGA11) which is induced across multiple models of reversible quiescence in distinct cell types and that has low expression in proliferating and senescent cells. Although ITGA11 is dispensable for proliferation of asynchronously cycling cells, it is required for efficient cell-cycle re-entry from quiescence. Mechanistically, ITGA11 promotes YAP accumulation and nuclear localization, thereby sustaining SKP2 expression and p27 degradation during cell cycle re-entry. Depletion of p27, or pharmacological activation of YAP signalling rescues the cell-cycle re-entry defect caused by ITGA11 depletion. Together, these findings identify ITGA11 as a functional quiescence signature protein that couples extracellular matrix sensing to YAP-dependent regulation of the Skp2-p27 axis, revealing a mechanism that controls the transition from quiescence to proliferation.

Abstract Imaging-based spatially resolved transcriptomics (imSRT) technologies provide high-throughput molecular-resolution spatial characterization of genes within cells. Conventional analysis methods to identify cell-types and states in imSRT data rely on gene count matrices derived from tallying the number of mRNA molecules detected for each gene per segmented cell, thereby overlooking subcellular heterogeneity that can be useful in defining cell states. To take advantage of the molecular-resolution information in imSRT data and potentially identify cell-states based on subcellular heterogeneity, we developed STARIT (Spatial Transcriptomics As Rasterized Image Tensors). STARIT converts transcripts within segmented cells in imSRT data into an image-based tensor representation that can be combined with deep learning computer vision models for downstream analysis. Using simulated and real imSRT data, we demonstrate that STARIT distinguishes transcriptionally distinct cell-types and further separates cell states based on subcellular transcript localization, which conventional gene count analysis fails to capture. By providing a standardized framework to encode subcellular molecular information in imSRT data, STARIT will enable deeper insights into subcellular heterogeneity and enhance the identification and characterization of cell-types and states that are overlooked by gene count representations.

Protein synthesis supports most biological processes. In the brain in particular, protein synthesis plays a critical role in physiological and pathological states. Here, we describe Tellurophene-Alkyne Cycloaddition-mediated Amino acid Tagging (TeACAT), a versatile strategy for fast, facile, and flexible tagging of newly synthesized proteins in mice. TeACAT is based on metabolic incorporation of the non-canonical amino acid TePhe into proteins by the endogenous protein synthesis machinery. Due to their high similarity, TePhe can efficiently replace canonical Phe without dietary or genetic manipulation. The subsequent bio-orthogonal reaction of TePhe with either fluorescent dyes or affinity handles enables both visualization and affinity enrichment of proteins synthesized during TePhe exposure. TeACAT is compatible with immunofluorescence for cell-type specific visualization of protein synthesis with subcellular resolution and can be used in conjunction with routine proteomics to identify and quantify newly synthesized proteins. Robust incorporation into the mouse proteome was observed on the scale of hours to days, allowing the interrogation of various biological processes. In summary, TeACAT enables the visualization and quantification of protein synthesis with minimal perturbation for biological discoveries.

Background: Fibrosis, characterised by excessive accumulation of collagen type I (COL1), is a common feature of chronic diseases, including liver diseases (LDs), chronic kidney disease (CKD) and heart failure (HF). COL1 degradation products can be detected in urine by proteomics/ peptidomics analyses and may serve as non-invasive biomarkers of fibrosis. We aimed to identify a common molecular signature of fibrosis across these diseases that may ultimately guide interventions to slow disease progression and prevent organ damage. Methods: Using capillary electrophoresis coupled to mass spectrometry (CE-MS), naturally occurring COL1 degradation products (peptides) in the urine of patients with fibrotic disease, LDs (n=127), CKD (n=263) or HF (n=187), were investigated and compared with the same number of matched controls. Disease-associated COL1 peptides were identified separately for each condition, and peptides showing consistent associations across the three diseases were selected to define a common fibrosis signature. A support vector machine model based on the selected peptides was developed and validated in independent cohorts of patients with LDs (n=110), CKD (n=93), HF (n=32) and controls (n=643). Results: We identified a common fibrotic signature consisting of 50 COL1 degradation products, mainly downregulated in fibrosis. A model based on these peptides achieved a strong performance, with an area under the receiver operating characteristic curve (AUC) of 0.935 (95% confidence interval (CI) 0.917-0.953, p<0.0001) in an external validation cohort comprising pooled disease groups (LDs, CKD, and HF) and controls. Performance was maintained in LDs, CKD and HF, with AUCs of 0.917 (95% CI 0.890-0.944, p<0.0001), 0.951 (95% CI 0.931-0.971, p<0.0001) and 0.950 (95% CI 0.903-0.997, p<0.0001), respectively. The model scores were significantly associated with fibrosis stage in LDs (p=0.0097) and with interstitial fibrosis and tubular atrophy in CKD (p=0.045). Conclusion: A model of urinary COL1 peptides captures a shared collagen degradation signature across organs and diseases, enabling the non-invasive assessment of fibrosis irrespective of its origin. As these peptides exclusively reflect collagen degradation, the findings suggest impaired collagen degradation as a driver in fibrosis. Future clinical studies are warranted to evaluate the utility of this model for early fibrosis detection and earlier implementation of anti-fibrotic interventions.

Metabolic dysfunction-associated steatotic liver disease (MASLD), the leading cause of chronic liver pathologies worldwide, represents a growing clinical burden. Its diagnosis remains reliant on liver biopsy that limits early detection and the ability to capture molecular changes across disease progression. A systematic understanding of stage-dependent gene expression changes is essential to identify biomarkers and effectively characterize disease mechanisms. Therefore recent studies provided databases for searching genes as well as prediction of multi-gene signatures for disease progression. However, there is still a need for interactive and comprehensive meta-analysis of datasets of MASLD patients with available histological metadata. Herein, we performed a meta-analysis of RNA-seq datasets using NAFLD Activity Score (NAS; n = 897) and fibrosis stage (n = 856) upon conducting pairwise comparisons across histological stages and identified differentially expressed genes associated with disease progression. Most importantly, we provide our findings via a dedicated web server, the MASLD-META NETWORK (https://masld.scilicium.com), enabling users to interactively explore meta-analysis results across diverse network modalities. In addition, we characterized gene expression dynamics across increasing disease stages to identify consistent progression-associated pathways using Louvain clustering. Network-based parameters such as centrality in combination with meta-analysis scores further highlighted central genes and pathways implicated in disease mechanisms. Accordingly, MASLD-META NETWORK enabled an integrative reassessment of recently published gene signatures, identifying COL1A1, COL3A1, THBS2, FBLN5, and PDGFA as the most central genes, and SULF2, MMP14, IL32, GPNMB, and COL3A1 as candidate markers of earlier transcriptional alterations. Network analysis of MASLD associated biological modules further identified LAMA2 and LAMA3 as previously unrecognized central candidate targets.

Lysine acylations, including lactylation (Klac), are pivotal regulators of cellular physiology. However, their analysis is currently bottlenecked by antibody enrichment strategies that suffer from sequence bias and require milligram-scale protein inputs, severely precluding the profiling of scarce clinical biopsies and rare cell populations. Here we present ChemIntelligence, an acyl-NHS chemistry-empowered derivatization strategy that rapidly generates unprecedented acylation-specific spectral libraries, exemplified by over 2.5x10^9 human Klac peptides, enabling cross-species reference atlases. Integrated with Prosit-based rescoring, these libraries substantially increase Klac identifications across diverse proteomic datasets. Leveraging this spectral resource, we devised ChemIntelligence Scope, a reproducible, multiplexed parallel reaction monitoring (PRM) platform that quantifies hundreds of Klac peptides per injection from as little as ~200 ng of cell lysates, clinical biopsies, and even true single cells - revealing functional Klac signatures inaccessible to conventional methods. The ChemIntelligence pipeline also extends seamlessly to lysine nicotinylation, underscoring its broad adaptability for discovering and profiling new acylations. Together, these chemical and computational advances establish a scalable, antibody-free framework for acyl-proteome mapping that overcomes input constraints and enables deep functional insights from otherwise intractable biological samples.

The Gram-negative bacterial cell envelope features an asymmetric outer membrane, that confers intrinsic resistance to toxins. Maintenance of this barrier relies on the Mla system, which mediates retrograde transport of mislocalized phospholipids. In Escherichia coli, this system comprises the lipoprotein MlaA, the periplasmic shuttle protein MlaC, and the ABC transporter complex MlaFEDB. Intriguingly, in Campylobacter jejuni, mlaA and mlaC share an operon with an encoded Resistance-Nodulation-cell Division antiporter potentially involved in anterograde phospholipid transport. Here, we describe the functional and mechanistical characterization of Cj MlaC. Complementation experiments in E. coli show that Cj MlaC functions independently of the native Mla system. Native mass spectrometry revealed that Cj MlaC uniquely exists as both monomer and dimer. Lipid binding stabilized the dimer and ion mobility mass spectrometry showed that conformational transitions precede phospholipid release, suggesting a cycle between a low-affinity monomer and a higher-lipid-affinity dimer. Cj MlaC binds phospholipid species distinct from Ec MlaC, showing an increased propensity for lysophospholipids, consistent with the unusually lysophospholipid-rich lipidome of C. jejuni, indicative of evolutionary adaptation to this unique lipid environment. Collectively, these findings uncover structural and mechanistic features of Cj MlaC and support divergent physiological roles for Cj and Ec MlaC in phospholipid trafficking.

Decellularized extracellular matrices (ECMs) are widely used in regenerative medicine, yet current evaluation criteria prioritize cellular removal rather than preservation of the ECM characteristics that govern tissue behavior. Here, we demonstrate that efficient decellularization is achieved across a broad range of chemical conditions, whereas preservation of structurally and biologically relevant ECM components is confined to narrow, tissue-specific windows defined by coupled detergent interactions. Quantitative proteomics revealed that intrinsic ECM composition is strongly associated with tissue-specific susceptibility to decellularization-induced damage and provided molecular context for the distinct preservation responses between tissues. Optimized matrices retained major structural ECM components and supported tissue-specific cellular organization and cell-mediated mechanical reinforcement following cellular repopulation despite uniformly low residual DNA across protocols. Together, these findings support a shift in decellularization quality assessment from DNA-based evaluation toward preservation of biologically relevant ECM and establish a composition-driven strategy for the rational design of regenerative biomaterials with tissue-relevant biological and mechanical properties.

Pseudomonas aeruginosa is an opportunistic human pathogen capable of infecting a wide range of tissues and organs. Its persistence during chronic infection is strongly associated with biofilm formation, which depends on extracellular polysaccharides such as alginate. The HtrA-like periplasmic serine protease MucD is a key regulator of bacterial virulence, stress response, and alginate production, yet its molecular mechanism has remained largely unclear. Here, we discovered the alginate acetylation and export proteins AlgX and AlgK as MucD substrates, and characterized their degradation by mass spectrometry and bioinformatic analysis. We further determined the cryo-EM structure of MucD bound to an AlgK-derived substrate peptide, offering atomic insights into MucD oligomerization assembly, substrate recognition, and specificity. Together with structure-guided mutagenesis and biochemical assays, our results revealed that MucD proteolytic activity is governed by an equilibrium between a resting 12-mer and an active trimer. Crucially, we demonstrate that MucD represses alginate biosynthesis post-translationally, in addition to its previously implicated role in transcriptional regulation. These findings define a distinct activation mechanism and regulatory function for MucD and provide new insight into bacterial HtrA-like serine proteases.

The immune checkpoint protein PD-L1 plays a pivotal role in tumor immune evasion by binding to PD-1 on immune cells, including T lymphocytes. While the expression and function of PD-L1 have been well studied, the importance of its spatial organization on the cell surface of tumor cells remains poorly understood. In this study, we used super-resolution microscopy combined with biochemical perturbations to investigate the factors regulating PD-L1 clustering and its effects on PD-1 binding and T cell inhibition. We found that PD-L1 is organized into nanoscale clusters at the plasma membrane, with distinct regulatory roles for the actin cytoskeleton, galectin-3, and cholesterol. Disruption of cortical actin increased PD-L1 cluster size, while galectin-3 promoted smaller, denser clusters and increased PD-L1 lateral mobility. Cholesterol depletion reduced PD-L1 cluster size and number and impaired PD-1 binding. These findings indicate that PD-L1 surface organization is collectively regulated by the actin cytoskeleton, galectin-3, and membrane cholesterol within the plasma membrane of tumour cells. Our results provide new insights into the dynamic regulation of PD-L1 and its potential as a therapeutic target in cancer immunotherapy.

Causal mediation analysis is widely used to identify biological pathways linking exposures to outcomes, but most methods assume homogeneous mediation effects across individuals. In high-dimensional omics settings, this assumption can mask important heterogeneity driven by demographic, genetic, or environmental factors. We propose the M-high-learner, a flexible framework for detecting heterogeneous mediation effects with high-dimensional mediators. The method identifies mediators with subgroup-specific indirect effects while distinguishing them from null or homogeneous signals and controlling the type I error rate. It is computationally efficient, scalable, and yields interpretable sub-types. Simulation studies show that the proposed approach achieves high power while maintaining accurate error control. Applications to the Framingham Heart Study and the Multi-Ethnic Study of Atherosclerosis reveal that the mediation role of gene expression in sexs effect on high-density lipoprotein varies across subgroups defined by body mass index and age. Our framework provides a practical tool for uncovering heterogeneous biological mechanisms in high-dimensional genomic studies. Author SummaryBiological processes linking risk factors to disease often differ across individuals, but many existing methods assume these processes are the same for everyone. This can hide important differences between groups. We developed a powerful method to identify when these pathways vary across subgroups using large-scale molecular data. Our approach detects differences in how intermediate biological factors contribute to outcomes in populations defined by characteristics such as age and body mass index. Applying our method to population studies, we found that some biological pathways operate differently across groups, suggesting that key mechanisms may be missed when differences are ignored. Our work provides a tool to better understand how disease-related processes vary across individuals, which may support more targeted and personalized approaches to health research.

Adhesion of cells is the key factor determining functioning of multicellular organisms. Viscoelastic properties of cells can be studied by multiple methods. However, attractiveness of cells or extracellular matrix without the elastic component (dispersive adhesion) is not accessible. We present an extension of force spectrometry technology: the Multivalent Adhesive Probe Atomic Force Microscopy (MAPA) that delivers dispersive adhesion maps of live cells and biosurfaces, and identifies differences unresolved by viscoelastic probing.

Genome-wide coexpression analysis of intact tissue samples is a powerful approach for identifying reproducible signatures of cell types and states, since it can survey vast numbers of individuals, cells, and transcripts. However, it can be difficult to optimize gene coexpression network construction and compare results from independent analyses. To address these challenges, we developed OMICON (theomicon.ucsf.edu) for research on human brain gene coexpression networks. OMICON contains gene expression data from >17K normal and neoplastic human brain samples with standardized metadata. Systematic analysis of independent datasets identified >250K gene coexpression modules, which were characterized and compared via enrichment analysis with >40K gene sets. All modules are discoverable via an advanced search engine that can filter by genes, metadata, and enrichment results. Analyses can also be browsed with an interactive workflow visualization tool, and users can communicate within OMICON using @mention functionality to support communal research on human brain gene coexpression networks.

Mutations in the E3 Ubiquitin (Ub) ligase RNF216 cause Gordon Holmes syndrome (GHS), a neurodegenerative disorder accompanied by neuroendocrine disruption. We developed an orthogonal ubiquitin transfer (OUT) platform to capture RNF216 substrates in neuronal cells and identified OTUD4, a deubiquitinating enzyme (DUB) mutated in GHS, and FMRP, a neuronal-enriched translational repressor. RNF216 predominantly synthesizes K6-linked Ub chains on OTUD4 to induce its degradation, forming donut-shaped structures in neurons. In return, OTUD4 removes the ubiquitination of RNF216 and FMRP. Analysis of RNF216 substrates revealed biological functions regulating protein synthesis, a shared function of the OTUD4-RNF216 substrate interaction network. Indeed, RNF216 expression increased protein synthesis rates in different cell types while Rnf216 deletion decreased dendritic development in neurons. Overall, our findings show that RNF216 and OTUD4 balance rates of protein synthesis and degradation and suggest GHS-related mutations in RNF216 or OTUD4 may offset this balance, triggering neurodegeneration.

Three-dimensional (3D) mural cell morphology is heterogeneous and coupled to vessel geometry, however, measurements from two-dimensional (2D) maximum intensity projections (MIP) obscure overlapping processes and cell-vessel contacts. Accordingly, we developed Mural-VISTA, a semi-automated Python workflow for mural cell-vessel interaction and single-cell topo-morphology analysis of reconstructed surface meshes. This workflow integrates mesh pretreatment, interactive centerline extraction, hierarchical segmentation of cell soma, main axis and secondary processes (branches), and extraction of 36 multiscale (cell process segment level, process level, and whole cell level) topo-morphological and vessel-referenced metrics. Mural-VISTA identified morphological changes in pericytes and vascular smooth muscle cells (vSMCs) with altered RhoA activity. Constitutive active RhoA (RhoA CA) over-expression reduced branch complexity and increased process alignment in both cell types, while increased whole-cell and branch solidity only in vSMCs. Dominant negative RhoA (RhoA DN) over-expression increased branch abundance and reduced branch solidity in pericytes but not vSMCs, suggesting cell-type specific effect of reduced RhoA activity. In conclusion, Mural-VISTA enables quantitative 3D profiling of mural cell architecture and its spatial relationship with the vessel.